using System;
using System.IO;

using Atomic.Libraries.Plotting;
using Atomic.Libraries.Plotting.PlanePlotting;
using Atomic.Vasp.LinearResponse;

public class Program
{
	public static void Main(string[] args)
	{
		// Parse VASP linear response output.
		VaspLinearResponseResult r = new VaspLinearResponseResult(new DirectoryInfo("/home/mbjba/results/Fe4N-20131003/Fe4N"));

		// Compute phonon density of states and plot.
		r.Phonons.StateDensity.Plot();

		// Compute phonon free energy and plot.
		Gnuplot.Plot(new AxisLabel("Temperature, K"), new AxisLabel("Energy per atom, eV"),
			new FunctionPlot(r.Phonons.ThermalPropertiesPerAtom.FreeEnergyFunction, 0.0, 2000.0, new PlotTitle("F_ph")));
	}
}
